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Abstract: The paper presents an electronic self-oscillation acoustic system (SAS) destined to measure of stresses 

variations in the elastic media. The system consists of piezoelectric detector, amplifier-limiter, pass-band 
filter, piezoelectric exciter and the frequency meter. The mass rock plays a role of delaying element, in 
which variations in stresses causing the variations of acoustic wave velocity of propagation, and successive 
variation in frequency of oscillations generated by system. 

The laboratory test permitted to estimate variations in frequency caused by variations in stresses of 
elastic medium. The principles of selection of frequency and other parameters of the electronic system in 
application to stresses measurement in condition of the mine were presented. 
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1. Introduction 

Autodyne systems are widely applied in radio engineering for monitoring of changes in tested object [5,6]. In 
mechanics and geophysics autodyne systems are used much rarely. Given paper describes a self-exiting acoustic 
system of autodyne type destined to measurements of stressed changes in elastic media. 

The structure of mentioned systems consists of amplifier with feedback loop of such properties, that the entirely 
system is unstable and generates oscillations. The function of sensitive element is performed by elements of 
feedback loop, which properties are depending on variations of measured value. As the final effect the frequency 
of self-oscillations is a measure of variations in input value. 
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2. Basic equations of SAS 

Lets consider self-oscillation acoustic system, presented in Fig. 1, which consists of 

• piezoelectric detector PD 

• amplifier-limiter A-L 

• pass-band filter F 

• piezoelectric exciter PE 

• delaying element DE 

• frequency meter FM 
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Figure 1. Basic schema of self-oscillating acoustic system which varies frequency under influence of stress in delaying element 



Input of A-L is a harmonic signal described by amplitude U, frequency uj and phase angles (pu'. 



u(t) = U cos{ujt + (f>u) 

Output signal v{t) of non- linear A-L contains the basic frequency as well higher harmonics: 



(1) 



v(t) = Fu;{v) cos(cjt -\- (fi) -\- F2uj{v) cos{2(jjt + (p2) + •.. 



(2) 



Amplitudes of harmonic components are depending on amplitude of input signal U. The most important is the 
first component of sum (3): 



v{t) = Fuj{v) cos{(jjt + (fi) 



(3) 



Non- linear characteristic of amplifier-limiter is shown in Fig. 2, under assumption, that an amplifier is saturated 
after increasing of amplitude U. 
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Figure 2. Non-linear dependence between output amplitude of A-L and input amplitude U 

Due to elimination of higher harmonic, signal v is transmitted by pass-band filter F, realised as parallel connection 
of inductivity L, capacity C and conductance G , which is shown in Fig. 3 
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Figure 3. Filtration of higher harmonics in the resonant circuit 



Admittance of such filter for the frequency ou can be written as: 



Y(uj) = iujC + —- + G, 
lujL 



(4) 



where the real part of the admittance G contains the inner conductance of current source Gz^ component Of 
describes energy loss in the filter, and finally Gw component describes energy loss in piezoelectric exciter: 



G = Gz + Gf + Gu 



(5) 



The output signal of A-L was presented as parallel connection of current sources, corresponding with successive 
harmonics cj, 2c(;, ...muj. The basic component of W{t) on the filter output can be described as 



W(t) = Wuj cos [Ujt + Lpn,{00)\ 



(6) 
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where denotes the absolute value, and (pw - phase angle of the complex amplitude: 

= W^e'""- (7) 
which is connected with amplitude of first harmonic Fi{U) by relationship: 

where Y(uj) denotes filter admittance. It implies that: 

~ \Y{,j)\ 

as well as : 

= - arg [Y{uj)] = (fu - arctan (^~~~q^^ (10) 
The complex amplitudes of higher harmonics are described using formulae analogous to (8) 

W^. = ^4^, (11) 

The absolute value of admittance Y{muj) is significantly greater than absolute value of admittance Y(uj) for values 

of cj, nearly to resonant frequency cjq = 1/V LC, for which Y = (cj) ~ G, so that : |Wma;| ^ |^cj|, which permits 

to omit the influence of higher harmonics, and to perform the analysis using only component (6). 

Voltage signal (6) drives the input of PE, exciting acoustic waves in delaying element DE, forming the feedback 

voltage: 

u{t) — kw{t — r) = kWco cos(cjt — cjr + ipw) (12) 

where: 

r - delay time, 

k-coefficient of damping and energy lossess in delaying element. 
Substituting (11) to (12) we obtain: 

u(t) = cos(cjt - ijOt + - arg Y{ijo)) (13) 

\^ \^)\ 

Comparing (13) with the formula describing input signal of the amplifier (1), we can to obtain the equation of 
phase balance: 

ujT + arg Y{uj) = 2iin (14) 
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which determines the frequency of self-oscillations. Simultaneously we have the following equation of amplitude 
balance, 

which permits to estimate the necessary gain factor of A-L element and to determine the amplitude of self- 
oscillations. 

3. Relationship between frequency of self-oscillation and 

stresses 

The dependence of left side of phase balance equation (14) on the frequency is shown in Fig. 4b, where: 



(p{(jo) = argy(c<;) = arctan 



iu;C-^)/G 

UjIj 



(16) 



The pass band of resonant circuit of Q factor Acj = {ujq ~ Iq'^o + Iq) cover the some frequencies ujn 
corresponding with the different values of index n. To avoid generation of multiple frequencies, the following 
inequality has to be fulfilled: 



Using the estimation 



CJn+l - OJn > (17) 



Un ^ (18) 

T 



which corresponds the omitting of component in (14) we have: 

1 > ^ Q > c^oT (19) 
r Q 

In such conditions frequency can be determined as follows: 

LOT = 2nn. (20) 
Deriving (20) with stress a as the independent variable we obtain: 

du dr ^ 1 duj 1 dr 

--^oj-— = 0^ ---r= 21 
dcr dcr cjdcr r da 

Formula (21) is a basic to determine the stresses variation using SAS. Lets take in account phase component in 
equation (14) by approximation: 
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, . d(^(cJo) ( X 2(5 \ -rr( \ 

^ — [uj - uo) — [uj - uJo) = T[yj - ujq) 



Where T = ^ 

UJQ 



^ denotes time constant of resonant circuit. 
Substituting (21) into (13) we have: 



(22) 



a)T + T(w - Wo) = 27rn 



(23) 



and: 



dcj, dr 
-—{r + T)+uj--^{) 
dcr do- 



As result, instead equation (21) we have equation: 



(24) 



1 duj 



1 dr 



(25) 



T + T do- 

where ^ component was changed into Usually T <^ r and for this reason derivative ^ in (25) is almost 

equal to (21). 



Janusz Kwasniewski, Yury Kravtsov, Ireneusz Dominik Lech Dorobczynski 



Equation (23) makes possible to estimate the sensitivity of variation in frequency caused by variation in temper- 
ature 



which allows to write: 



T 



As result we obtained that, the sensitivity of variation in frequency caused by variation in temperature is 
times lower than the sensitivity of variation in resonant frequency caused by variation in temperature. 

4. Estimation of SAS sensitivity 

Computer simulations performed in MATLAB-Simulink environment were based on preliminary research 
performed in laboratories of AGH University of Science and Technology in Krakow. The mentioned research was 
concerned on determination of dependency between compressing stress and velocity of sound in the specimen 
made of concrete. Obtained linear regressive model has a form: v = a ■ F -\- b 



where: 

F - compressing force expressed in kN, 
V - sound velocity expressed in m/s. 

Values of coefficients a and b are as follows: a = 5.27807339449537 [m/Ns], b = 2678.4917 
Obtained results are presented in Fig 5. 
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Figure 5. Dependency between velocity of sound and pressing force in concrete 

Simulation was performed for specimen length L=0.45 m. Parameters of model built in Simulink were as below: 
• static characteristic of limiter was described by formula y = arctan(x). 



7 



Self-oscillation acoustic system destined to measurement of stresses in mass rocks 



gain factor of amplifier was assumed as K—25 

band pass filter, consisted of two sections, each of them described by transfer function: 



Gf{s) 



(27) 



where: 

= 30787.61 ^ 



UJ2 



— 30787.61^ - are lower and upper cut-off frequencies of the filter. Obtained resonant 



frequencies in dependence of pressing force are presented in Table 4 



F [kN] 


9.81 


19.62 


29.43 


39.42 


49.05 


V [m/s] 


2717 


2768 


2819 


2871 


2922 


/o [Hz] 


4666 


4676 


4688 


4694 


4702 



Basing on previous data, system sensitivity factor, expressed in form: 



can be estimated as: 



S= 1.02 



Hz 
kN 



(28) 



5. Laboratory stand 



During our preliminary research at the Department of Process Control at the AGH University of Science and 
Technology the Self-exited Acoustical System has been developed. System diagram is shown in figure 6, where 
amplifier, shaker (E) and receiver (R) are formed in a feedback loop. The shaker (E) was fixed to a stone bar. 
The bar was put into a testing machine to create a load. On the beams surface four accelerometers were fixed: 
three of them on the same surface as the shaker and the last one directly on the opposite site. 
As a result of positive feedback, there is a bilateral interaction between the control device and the vibrating 
system, which allows the self-oscillating system to control its own energy balance. As a result of that, despite the 
fact of losses of energy occurring in the system, there are non-fading periodic oscillations. Many applications of 
the self-oscillating phenomenon are known. These are: the vibrations of cutting tools, turbine blade vibration, 
vibration of aircraft wings, vibration of bridge suspension, which may cause a destruction of the bridge, such as 
the Tacoma Bridge, etc. In these systems, we are trying to eliminate these vibrations. The system was intended 
to measure a change of stress in elastic mechanical structures, construction and rock masses. The purpose of 



Janusz Kwasniewski, Yury Kravtsov, Ireneusz Dominik Lech Dorobczynski 




Figure 6. Self-exited Acoustical System diagram, where E - shaker and R - receiver 

this study is to determine the possibihties of using this system for real objects such as bridges, dams, buildings, 
mines, etc. The sensitivity of this system, for small and large deformation, is higher than the sensitivity of 
other measurement systems (especially open systems). Tests on a single sample of sandstone were performed. 
They tried to examine the impact on the stress measurement parameters such as: position of sensors, position of 
actuator, and the influence of geometrical shape and dimensions of sample. 




Figure T. The test stand for the endurance tests using the SAS System 

The study used a sample of sandstone compressed in the frame by a hydraulic press (Fig. 7). In addition, the 
system was equipped with a force sensor to calculate the compressive stress in the beam in real time. Accelerations 
were measured by three accelerometers for three reasons: 

- to determine what impact the distance from the emitter to the receivers has on the results and the position of 

resonances, 

- to avoid a situation where the sensor is located in a resonance node, which would result in incorrect test results, 

- to designate the velocity of wave propagation in a material with the correlation methods. 
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The study was conducted in three configurations of the position of emitter (E). For a sandstone sample with 
rectangular cross-section of 60x70 emitting and receiving devices were mounted in the configuration shown in 
Figure 8. 
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Figure 8. Diagram of deployment of emitters on sandstone 



6. Experimental results 

The SAS is designed for measurements of stress in elastic structures including rocks as well as mines, bridges, 
dams, buildings, etc. Therefore it is necessary to create a universal procedure to detect the change in the strain 
in objects mentioned above. The designers decided to use first the system to conduct tests in the open loop 
[3]. To the sample, charged with various compressive forces, was initially given the chirp signal emitted in the 
frequency range of 100 Hz to 20 kHz. This allowed to obtain the amplitude - frequency characteristics (Fig. 9) of 
the structure with visible resonances (Fig. 10). 

Figure 10 presents the movement of the first resonance peak believed to be related to the change of stress in the 
sample. Changing the frequency respectively equals to 90 Hz. For the closed loop system, which is shown in figure 
11, the frequency change at the same load change is as large as 450 Hz. In addition, the amplitude of vibration 
is 20-30 times greater than in an open loop system, which makes them more visible and easier to identify. 
Thus our experiments have demonstrated an opportunity to detect stress changes in elastic geophysical objects: 
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Figure 9. Vibration transmissibility functions for different sample loads 
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Figure 10. Selected and expanded resonance peaks related to the change of stress 



marble, sandstone and also concrete by means self-oscillating acoustic system (SAS). We hope that similar system 
might be helpful for stress change detection in rock mass and mines as well as civil engineering objects i.e. bridges, 
dams and buildings. 
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Figure 11. Selected and expanded resonance peaks of closed loop system related to the change of stress 



6.1. Conclusions 



The main thesis presented in this paper is that self-exited acoustical system with the positive feedback allows 
monitoring of the change of the stresses in elastic structures. The usefulness of this phenomenon for different 
states of strain and the various sizes of test samples of sandstone, marble and concrete was examined to prove this 
property. This study confirmed the existence of the self-excitation phenomenon, which can be used to recognize 
the state of rock mass changeability. 
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